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Abstract 

In this paper, following |T1 [21 SI O E] we consider the relations 
between well-known Fourier transform algorithms. 
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1 Introduction 

A semifast algorithm for a computation is an algorithm that significantly 
reduces the number of multiplications compared with the natural form of 
the computation, but does not reduce the number of additions. A semifast 
Fourier transform in GF(q) is a computational procedure for computing the 
rz-point Fourier transform in GF(q) that uses about nlogn multiplications 
in GF(q) and about n 2 additions in GF(q) [I]. 

The n-point Fourier transform over GF(2 m ) is defined by 
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where W = i, j G [0,n — 1], is a Vandermonde matrix, 

an element a G GF(2 m ) of order ord(a) = n \ 2 m — 1 is the Fourier transform 

kernel. 

Consider a few well-known Fourier transform algorithms [U [21 El IH El El Ej 
on examples. 

The finite field GF(2 3 ) is defined by an element a, which is a root of 
the primitive polynomial x 3 + x + 1. Cyclotomic cosets modulo 7 over 
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GF(2) are (0), (1,2,4), (3,6,5). The binary conjugacy classes of GF{2 3 ) are 
(a ), (a 1 , a 2 , a 4 ), (a 3 , a 6 , a 5 ). The standard basis of GF(2 3 ) is (a , a 1 , a 2 ). 
The normal basis of GF{2 3 ) is {(3 1 , (3 2 , (3 4 ) , where (5 = a 3 . Another normal 
basis of GF{2 3 ) is the cyclic shift of the previous normal basis: (7 1 ,7 2 ,7 4 ), 
where 7 = a 6 . 

Consider the Fourier transform of length n = 7 over the field GF(2 3 ). Let 
us take the primitive element a as the kernel of the Fourier transform. The 
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Fourier transform of a polynomial f(x) = fix 1 consists of components 

Fi = /(«*), ie [0,6]. 
Thus, 
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2 Goertzel's algorithm (1958) 

Blahut's modification for finite fields (1983) 

The first step of the Goertzel algorithm is a long division f(x) by every 
minimal polynomial: 

f(x) = (x+l)q (x) + r 00 

f(x) = (x 3 + x+ l)gi(ar) + r^x), n(x) = r 21 x 2 + r u x + r m , 
f(x) = (x 3 + x 2 + l)q 2 {x) + r 2 {x), r 2 (x) = r 22 x 2 + r 12 x + r 02 

where 

^oo = fo 

roi = fo + fa + h + h 
rn =/i + / 3 + /4 + / 5 
r 2 i = f 2 + U + h + fe 
ro2 = fo + h + fa + h 
ri2 =fi + fa + fa + fa 
r 22 = fa + h + fa + fa- 
The second step of the Goertzel algorithm is an evaluated remainder in 
every element of the finite field: 
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r 2 iQ; 4 + r n a 2 + r i 



»"01 



= r 22 a 6 + r 12 a 3 + r 02 
= r 22 a 5 + ri 2 a 6 + r 02 
= r 22 a 3 + ri 2 a 5 + r 02 
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3 Blahut's algorithm (2008) 

We have 
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4 Fedorenko— Trifonov's algorithm (2002) 

A polynomial L(x) = ^^liX 2 \ k G GF(2 m ), is called a linearized poly- 

i 

nomial. Any polynomial can be decomposed into a sum of linearized poly- 
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nomials and a free term. A polynomial f(x) = fax 1 can be represented 



as 



i=0 



f(x) = L (x°) + L^x 1 ) + L 2 (x s ) 

L (y) = fo 

Li(y) = hy + hy 2 + hy 4 

Hy) = hy + hy 2 + hy 4 . 



We have 



L (a°) + L 1 (a°) + L 2 (a ) 
L (a°) + L^a 1 ) + L ' ' 



' 2 (« 3 ) = 



f(a 2 )= L (a°) + L 1 (a 2 ) + L 2 (a 6 ) = 
/(« 3 ) = 
/(« 4 ) = 



= L (l) + L 1 (a) + L 2 (l) + L 2 (a) 
L (l) + L 1 (« 2 ) + L 2 (l) + L 2 (a 2 ) 



/(« 5 ) = 
/(« 6 ) = 

Usinj 



L (a°) + ^(a 3 ) + L 2 (a 2 ) = L (l) + Li(l) + L^a) + L 2 (a 2 ) 
L («°) + Li(a 4 ) + L 2 (a 5 ) = L (l) + Ha) + Li (a 2 ) + 

L 2 (l) + L 2 (a) + L 2 (a 2 ) 
L (a°) + L^a 5 ) + L 2 (a l ) = L (l) + Li(l) + Li(a) + Li(a 2 ) + L 2 (a) 
L (a°) + Ha 6 ) + L 2 (a 4 ) = L (l) + Li(l) + L^a 2 ) + L 2 (a) + L 2 (« 2 ). 

Fi = f(a l ), these equations can be represented in a matrix form as 
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5 Trifonov— Fedorenko's algorithm (2003) 

A circulant matrix, or a circulant, is a matrix each row of which is obtained 
from the preceding row by a left (right) cyclic shift by one position. If entries 
of a circulant are matrices, the circulant is referred to as a block circulant. 
We call a circulant where the first row is a normal basis, a basis circulant. 

Let (P\ (3 2 , 13 A ) = (a 3 , a 6 , a 5 ) be a normal basis for GF(2 3 ). 

Using 
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6 Fedorenko's algorithm (1) (2006) 

Combining 
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we obtain the equivalent Fourier transform 



/F \ 




[ a 


a 


a 


a 


a 


a 


a 


F, 




a 


a 1 


a 2 


a 4 


a 3 


a 6 


a 5 


F 2 




a 


a 2 


a 4 


a 1 


a 6 


a 5 


a 3 


F± 




a 


a 4 


a 1 


a 2 


a 5 


a 3 


a 6 


Fs 




a 


a 3 


a 6 


a 5 


a 2 


a 4 


a 1 


F 6 




a 


a 6 


a 5 


a 3 


a 4 


a 1 


a 2 


W 




a 


a 5 


a 3 


a 6 


a 1 


a 2 


a 4 



fl 
h 
h 



h 
h 



1 


1 


1 


1 


1 


1 


1 " 




" 1 








ffo\ 


1 





1 


1 


1 












P 1 


P 2 


P 4 










fl 


1 


1 





1 





1 









P 2 


P A 


P 1 










h 


1 


1 


1 











1 






P A 


P 1 


P 2 










h 


1 


1 








1 





1 












P 1 


P 2 


P A 




h 


1 





1 





1 


1 















P 2 


P A 


P 1 




h 


1 








1 





1 


1 












P 4 


P l 


P 2 




w 



The equivalent Fourier transform F e = W e f e has the following structure: 



F e = WJ e = A e D e f e 
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where A e is a binary matrix, 

A e consists of binary circulants A^, 

D e is a block diagonal matrix, 

D e consists of basis circulant matrices Cj, 

/ is the number of cyclotomic cosets modulo n over GF(2). 
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7 Fedorenko's algorithm (2) (2006) 

Let (7 1 ,7 2 ,7 4 ) = (a 6 , a 5 , a 3 ) be a normal basis for GF(2 3 ). 
Combining 
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Note that the binary matrix A e for this algorithm in ([T]) consists of binary 
circulants and block circulants. 



8 Complexity 

The Fourier transform algorithms [SJ El E] of length n = 2 m — 1 over GF(2 m ) 
take two stages: 

1. The first stage is calculation of I m-point cyclic convolutions; 

2. The second stage is multiplying the binary matrix A e by the vector 
DJ e . 

The complexity of the first stage is about nlogn multiplications and 
additions over elements of GF(2 m ). The complexity of the second stage is 
N a dd < 2n 2 /logn additions over elements of GF(2 m ) [S]. 
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